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Abstract. We discuss the origin of the microscopic description of correlations in 
quantum many-particle systems obeying Fermi-Dirac and Bose-Einstein statistics. 
For correlation operators that give the alternative description of the quantum state 
evolution of Bose and Fermi particles, we deduce the von Neumann hierarchy of 
nonlinear equations and construct the solution of its initial-value problem in the 
corresponding spaces of sequences of trace class operators. The links of constructed 
solution both with the solution of the quantum BBGKY hierarchy and with the 
nonlinear BBGKY hierarchy for marginal correlation operators are discussed. The 
solutions of the Cauchy problems of these hierarchies are constructed for initial data 
satisfying a chaos property. 

Key words: von Neumann hierarchy; Fermi-Dirac and Bose-Einstein statistics; 
correlation operator; density matrix; marginal operators; quantum many-particle 
system. 



2000 Mathematics Subject Classification: 35Q40; 47D06. 



Dynamics of Correlations 



2 



Contents 

1 Introduction 

2 The von Neumann hierarchy for correlation operators of Fermi and Bose 
many-particle systems 3 

2.1 Preliminary facts [3 

2.2 Cluster expansions of density operators 

2.3 The von Neumann hierarchy for correlation operators 

2.4 A formal solution of the von Neumann hierarchy 

2.5 Correlation operators of particle clusters [13 

3 The Cauchy problem of the von Neumann hierarchy in the spaces of se- 
quences of trace-class operators |lj 

3.1 Some properties of a solution l4 

3.2 The group of nonlinear operators generated by the von Neumann hierarchy ... 16| 

3.3 The existence and uniqueness theorem lol 



4 On applications of the von Neumann hierarchy |21 

4.1 Functionals of the mean- value and the dispersion of observables 21 

4.2 Marginal density operators 24 

4.3 Marginal correlation operators 25 

4.4 On solutions of the BBGKY hierarchies in case of chaos initial data |26 

5 Conclusions |27 
References 28 



1 Introduction 

The paramount importance of the mathematical description of correlations in numerous prob- 
lems of the modern statistical mechanics is well-known. We refer, for example, to such fun- 
damental problems as quantum measurements and the rigorous derivation of quantum kinetic 
equations [IHS], in particular, kinetic equations that describe Bose gases in the condensate 
state [lOHll]. 

There are several approaches to the description of fluctuations in many-particle systems. 
Among them we mention the well-known approach of dynamics of correlations [TSIIIS] and its 
applications in plasma theory. Equilibrium states in framework of correlation operators are 
rigorously described in works [TTtlTS]. 

The aim of the paper is to formulate rigorously the evolution equations that describe non- 
equilibrium correlations of Fermi and Bose many-particle systems and to construct a solution of 
the corresponding Cauchy problem. Then using the constructed solution to establish its links 
with the solutions of the BBGKY hierarchies for marginal density and marginal correlation 
operators [T9H21]. 
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The paper is organized in the following way. 

In section 2 the description of evolution of states of Fermi and Bose many-particle systems in 
framework of correlation operators is given. The correlation operators are defined by means of 
cluster expansions of the density operators (density matrices), which are solutions of the initial- 
value problem of the von Neumann equations of a non-fixed, i.e. arbitrary but finite, number 
of particles. We deduce the von Neumann hierarchy that governs the evolution of correlation 
operators of Fermi and Bose particles with the general type of an interaction potential and 
construct a solution of its initial- value problem, in particular in case of initial data satisfying a 
chaos property. 

In section 3 we give the rigorous justification of the derivation of evolution equations for 
the correlation operators. We prove that the solution of the initial-value problem of the von 
Neumann hierarchy generates a group of nonlinear operators of the class Cq iii the suitable 
spaces of sequences of trace class operators. In these spaces the existence and uniqueness of a 
strong and a weak solution of the Cauchy problem of the von Neumann hierarchy is proved. 

In section 4 using introduced correlation operators we define functionals of average values 
and a dispersion of observables. Moreover, we give the description of states of Fermi and Bose 
many-particle systems in framework of marginal density and marginal correlation operators, 
that are adopted for the description of infinite-particle systems. The solutions of the correspond- 
ing BBGKY hierarchies for these operators are deduced on the basis of constructed solution 
of the von Neumann hierarchy. We demonstrate that the concept of cumulants of groups of 
operators of the von Neumann equations forms the basis for the solution expansions of men- 
tioned hierarchies of evolution equations. In the final section we give concluding comments and 
remarks. 

2 The von Neumann hierarchy for correlation operators 
of Fermi and Bose many-particle systems 

2.1 Preliminary facts 

Let y, be the one-particle Hilbert space, then the n-particle space Tin = H*^" is a tensor 
product of n Hilbert spaces "H. We adopt the usual convention that Ti®^ = C. Let (3„ be the 
permutation group of the set (1, . . . , n). We assign to each permutation tt G an isomorphism 
Pt^ of "H*^" onto itself. The operator p^^ transforms factorized elements ■0i ® '!/'2 ® • • • ® "'An ^ 
-^(g)n jj^^Q -0^(1) ® '07r(2) ® • • • ® '07r(n) ^ "H®". The symmetrizatiou operator and the anti- 
symmetrization operator S~ on "H®" are defined by the formula 



where |7r| denotes the number of transpositions in the permutation vr. The operators are 
orthogonal projectors, i.e. {S^y = S^, ranges of which are correspondingly the symmetric 
tensor product T-L^ and the antisymmetric tensor product of n Hilbert spaces T-L. We 
denote by J-"^ = 0^o^ri ^^e Fock spaces over Hilbert space 7/, that are associated with the 
corresponding Fermi and Bose systems [221123] . 




(1) 
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Hereinafter we consider a quantum system of a non-fixed, i.e. arbitrary but finite, number 
of identical (spinless) particles with unit mass m = 1 in the space W , v > 1, that obey Fermi- 
Dirac or Bose-Einstein statistics. The Hamiltonian H = system is a self-adjoint 

operator (Hq = 0) with domain V{H) = {ip = ^ipn e \ ipn e V{Hn) € Ht, En ll^n^nP < 
oo} C -F^. On functions ^„ that belong to the subspaces S^Ll{W) C V{Hn) C S^L'^{W) 
of infinitely differentiable functions with compact supports n-particle Hamiltonian Hn acts 
according to the formula 

n n n 

i^nV'n = J]i^W^n + 5^ J2 ^''"^H, ■ ■ ■ , tk)i^n, (2) 
1=1 k=l ji<...<ifc=l 

where K{i)il)n = —^Ag.ipn is the operator of the kinetic energy, 27ih is a Planck constant and 
. . . , ik)iJn = ^^''H<lh, ■ • • , <lik)^n is a fc-body interaction potential. We assume that the 
functions <l>(^\ k > 1, are translation-invariant, symmetric with respect to permutations of their 
arguments and satisfy the Kato conditions [2lj, which guarantee self-adjointness of operator 

(ED- 

States of systems of Bose and Fermi particles belong to the corresponding spaces = 
©J^O'^^('^n ) of sequences / = (/o, /i, ...,/„,. . .) of trace class operators /„ = /„(!, . . . , n) G 
2^(71^) and /o G C, that satisfy the symmetry condition: . . . , n) = fn{H, ■ ■ ■ jin) for 

arbitrary (^i, . . . , G (1, . . . , n), equipped with the norm 

oo oo 

"•^II^H-^^) = 5Z "-^"II^HW^) =Z1 Tri,...,„|/„,(l,...,n)|, 

n=0 n=0 

where Tri^...^„ are partial traces over 1, . . . , n particles [201123] . We denote by 2,1 the everywhere 
dense set in 2^{J-'^) of finite sequences of degenerate operators [23] with infinitely differentiable 
kernels with compact supports. Note that the space contains more general sequences 

of operators than those that determine states of systems. 

The Bose-Einstein and Fermi-Dirac statistics endows the states with additional symmetry 
properties. We illustrate them on kernels of operators [20|. Let . . . , Qn', Qi, ■ ■ ■ , Qn) is the 

kernel of the operator G 2\S^L'^{W)). In case of the Bose-Einstein statistics the kernel is 
a symmetric function with respect to permutations in each group of arguments 

fniQl, (l2, ■ ■ ■ , Qn', Ql, ^25 • • • > Qn) = fn{Qn{l) , (ln{2) , ■ ■ ■ , QTT{n)', Qn' (1) ^ ll' (2) ■> ■ ■ ■ 5 Qn'in))^ 

and in case of the Fermi-Dirac statistics the corresponding kernel is an antisymmetric function 
where tt, tt' G 

The evolution of all possible states, i.e. sequences D{t) = (1, Di{t, 1), . . . , Dnit, 1, . . . , n), 
. . .) G 2^{J-'^) of the density operators Dn{t),n > 1, is described by the initial-value problem 
of a sequence of the von Neumann equations (the quantum Liouville equations) [19 ] [20 ] [25] 

j^D{t) = -UD{t\ (3) 
D(t)|i=o = I^(0), (4) 
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where for / E ^l{J^u) ^ ^(-^) ^ ^\J^n) ^^e generator A/" = ®n=o^n of the von Neumann 
equation is defined by 

(Ar/)„ = --^ =-l(/„if„ (5) 

where Hn is the Hamihonian ([2]). 

In the spaces of sequences of trace class operators for initial- value problem (l3])-(j4]) 

the following statement is true. 

Proposition 1. A unique solution of initial-value problem is determined by the formula 

D{t)=g{-t)D{0), (6) 

where G{-t) = (B^^oGni-t), 

^n(-t)in = e-i*^"/.e^*^". (7) 

For D{0) e S^liJ^^) C 

is a strong (classical) solution and for arbitrary D{0) G 
2^{J-'^) it is a weak (generalized) solution. 

For the case of the Maxwell-Boltzmann statistics an analogue of this proposition is proved 
in [22], [23]. 

In the spaces 2.^{J-'^) the mapping : t — )■ G{—t)f is an isometric strongly continuous group 
that preserves positivity and self-adjointness of operators. 

For / e £o(-^k) T){—J\f) in the sense of the norm convergence in the spaces £,^{J-'^) there 
exists the limit by which the infinitesimal generator of the group of evolution operators ([7]) is 
determined: 

lim ^(^(-t)/ - /) = -^{Hf - fH) = -AT/, (8) 

where H = ©^q-^™ Hamiltonian ([2]) and the operator {—j^{Hf — fH)) is defined in the 

domain T>{H) d J-'^. 

Symmetrization and anti-symmetrization operators ([T]) are integrals of motion of the von 
Neumann equation and as a consequence 

and 

KS^ = s^K, 

i.e. the symmetry of states is preserved within evolutionary process. 
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2.2 Cluster expansions of density operators 

We introduce operators g{t) = (0, gi{t,l), . . . , gs{t,l, . . . , s), . . .) G 2}{J-'^) that are defined by 
the cluster expansions of the density operators D(t) = (1, Di{t, 1), . . . , Ds{t, 1, . . . , s), . . .) G 

Ds{t,Y) = gs{t,Y) + Yl l[g^x^^{t,X,), s > 1, (9) 

P:Y = [J^X„ X,CP 

|P| > 1 

where J2p y=[J- Xi \p\>i possible partitions P of the set y = (1, . . . , s) into 

|P| > 1 nonempty mutually disjoint subsets Xj C Y. 
We give instances of recursion relations (Q 

D,{t,l) = g,{t,l), 

D^it, 1, 2) = g^it, 1, 2) + Sig,{t, l)g,{t, 2). 

The sequence of operators g{t) defined by recursion relations ([9]) has a meaning of a sequence 
of correlation operators of Fermi and Bose many-particle systems, i.e. a sequence of operators 
that characterize the correlations of particle state^. We remark that cluster expansions ([9]) 
without operators Sf correspond to correlation operators that describe states of systems obey- 
ing the Maxwell-Boltzmann statistics. 

In order to construct the solution of recursion relations (Q, i.e. to express the correlation 
operators in terms of the density operators, we introduce the following notions. On sequences 
of operators /, / G 2^{J-'^) we define the ©-product 

(/ ® f)\Y\{Y) = Sfy^ f\z\{Z) J\Y\z\iY\Z), (10) 

ZCY 

where the operators iS|y| are defined by formula ([T]) and '^zcy subsets Z 

of the set Y = (1, . . . , s). We remark that similar product was introduced in [TTJIIH] for the 
studying of equilibrium states. 

By means of definition (ITU]) of the ©-product we introduce the mapping Exp^, and the inverse 
mapping Ln® on sequences h = (0, /ii(l), . . . , hn{l, . . . , n), . . .) of operators hn G S^iTi^) by 
the expansions 

(Exp^/^)|y|(r) = {1 + Y = (11) 

71=1 

P:y=Ui^> X.GP 

where 1 = (1, 0, . . . , 0, . . .) and 5|y|,o is a Kronecker symbol, and correspondingly, 

(Ln^(l + /i))|y|(y) = (5^(-l)"-i— )|^|(y)= (12) 

n=l 

= E (-i)"''"'(ipi-i)!'5^ n v.iTO- 

^■Y=yj,Xi x.dP 



The motivation of introduction of these operators is considered in section 4. 
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Then in terms of sequences of operators recursion relations ([9]) are rewritten in the form 

D{t)=Exp^ git). 

From this equation we obtain 

(/(t) =Ln®D(t), 

where D{t) = 1 + (0, D,{1), ^2(1, 2), . . . , Z^„(l, . . . , n), . . .)• 

Thus, according to definition ( |T0|) of the ©-product, in the component-wise form solutions 
of recursion relations ([9]) are represented by the expansions 

gs{t,Y) = D,{t,Y) + 

+ E (-i)"''''(ipi-i)!'5f n ^^1- (13) 

|p| > 1 

In terms of kernels of operators the first two equations of expansions (fT3|) have the following 
form 

9i{t,qi;q'i) = Di{t,qi;q[), 

g2{t, qi, q2] q'l, ^2) = ^2(i, qi, q2] q'l, ^2) - 

gi; q'i)Di{t, q2] ^2) ± Di{t, q2] q[)Di{t, qi, q'^)). 

The structure of expansions f|T3|) means that the correlation operators have a sense of cumulants 
(semi-invariants) of the density operators governed by the Cauchy problem of the von Neumann 
equation (^-@. 

Thus, correlation operators ( !T3|) give an alternative approach to the description of states of 
the Fermi and Bose many-particle systems, namely in the framework of dynamics of correlations. 

2.3 The von Neumann hierarchy for correlation operators 

From expansions f lT3|) using the von Neumann equation (|3]) we derive the hierarchy of evolution 
equations for the correlation operators g{t) = (0, gi{t, 1), . . . , gs{t, 1, . . . , s), . . .) G £,^{J-'^) 

j^g^itX) = ~K{Y)gs{t,Y) + (14) 

P:Y = \J^X„ ZiCXi, Z|p|CX|p|, XiCP 
|P|>1 ^17^0 2|P|7^0 

where 

Af^:^{l,...,n) = -'-[. ,^^-\l,...,n)], (15) 
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is a n-body interaction potential defined in ([2]), Yl p -y = \j-Xi,\p\>i ^^e sum over all 
possible partitions P of the set y = (1, . . . , s) into |P| > 1 nonempty mutually disjoint subsets 
Xi <zY and Ylz cx z-^Q ^ nonempty subsets Zj C Xj. We refer to the hierarchy 

of equations f lT^ as the von Neumann hierarchy. 

Hierarchy (|T4l) together with the initial state defined by 

9s{t,Y)l^^ = gMy), s>l, (16) 

describes the evolution of all possible states of Fermi and Bose many-particle systems in frame- 
work of correlation operators. 

We give an explanatory comment concerning the origin of the von Neumann hierarchy f|T^ . 
As a result of the differentiation of the definition of correlation operators (ITS]) by time variable 
in the sense of pointwise convergence in the spaces £,^(71^) and using cluster expansions (jH]) we 
obtain 

f^9s{t,Y)= (-1)'"'"'(|P|-1)!'5^ n ^1^.1(^.^0 = 

P:Y = \J^Xi XiCP XiCP 

= E (-l)"''"'(|P|-l)!'5^ E(-^^-l)™ n E U 9\z,At^Z,J. 

P:Y = \J^X, X,CP X,CP P' : X, = IJ^, ^j. ^J, CP' 

Changing the summation indexes P and P' in this expression and collecting similar terms gives 
us the chain of equalities ( IT^ which we interpret as the hierarchy of evolution equations for 
the sequence of correlation operator^. 

We cite an instance the typical equations of hierarchy f|T4|) 

^(?l(^,l) = -M(l)^/l(^,l), 

j^g2{t, 1, 2) = -Ar2(l, 2)g2{t, 1, 2) - ^^^^{1, 2)Sig,{t, l)g,{t, 2), 

j^g,{t, 1, 2, 3) = -^(1, 2, 3)g,{t, 1, 2, 3) - 

-{M^^h^,'^)+^^nh^,^)+^Sh^,'^,^))St9i{tA)g2{t,2,3)- 
-{Af^;;{l,2)+Af^;;{2,3)+AfShl,2,3))Stg,{t,2)g2{t,l,3)- 
-{ArS^{l,3)+Ar^;h2,3)+ArShl,2,3))Sig,{t,3)g2{t,l,2)- 
-Af^^h^, 2, 3)St{l, 2, 3)g,{t, l)g,{t, 2)g,{t, 3). 

In case of a two-body interaction potential the von Neumann hierarchy (fH|l reduces to the 
following form 

^^g,{t,Y) = -K{Y)gs{t,Y) + (17) 

+ E E E (-■^in?(^l'^2))5,Vi|(^'^lh^2|(^,^2), S>1, 

P:Y=XiiJX2 iiGXi 126X2 



^The rigorous justification is given in section 3. 
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where ^p.y=Xi u^2 possible partitions P of the set y = (1, . . . , s) into two 

nonempty mutually disjoint subsets Xi (Z Y and X2 C Y . 

We note that for classical many-particle systems hierarchy ( IT7|) was introduced in [26j as the 
zero-order approximation of the nonlinear BBGKY hierarchy for marginal correlation functions. 

In terms of kernels of correlation operators the first two equations of the von Neumann 
hierarchies f ll4p and fll7p have the form 

d 

d 

ifi-g^g2{t,qi,q2;q'i,q'2) = 

= ( - y E^^'^' - ^:) + - - '^^'^^'^'1 - q2)))92{t,quq2;q[,q'2) + 

1=1 

+ _ g^) _ ${2)(g^ _ ^^)) q^. q'^)g^{t, gs; q'2) ± 9i{t, q2] q'i)gi{t, qv, q'^)), 

where is a two-body interaction potential (|2]). 

We remark that a particular solution of the steady von Neumann hierarchy is a sequence of 
the Ursell operators, for example, its first two elements have the form 

9,(«,l) = e-"*-W, 

9,(U.2) = S±e-".S"'"(e-««''")-/), 

where /? is a parameter inversely proportional to temperature, the operators K(i) and are 
defined by We use the name Ursell operators to indicate their intimate structure with the 
expansions of classical equilibrium states developed by H.D. Ursell P?] . 



2.4 A formal solution of the von Neumann hierarchy 

To construct a solution of the Cauchy problem (fT4|) . ( [T6|) we first consider its structure for one 
physically motivated example of initial data, namely initial data satisfying a chaos property 
(statistically independent particles) [2T] . A chaos property means that there are no correlations 
in a system at the initial instant. In this case the sequence of initial correlation operators is 
the one-component sequence 

g{0) = {0,g,{0,l),0,...). (18) 

In fact, in terms of the sequence of the density operators this condition means that D{0) = 
(1, Di(0, 1), 52±Di(0, l)Di(0, 2), . . . Uti ^(0, 0, • • •)• 

On basis of representation f ll3p of the correlation operators in terms of the density operators 
and formula (E]), in view of the fact that the operators commute with groups of operators 
([7]), we have 

s 

9s{t,Y)= Yi (-i)"''"'(ipi-i)! n ^i^»i(-^'^^)'5^n^i(0'^)' ^^1- 

P:Y = lj^X, XiCP i=l 
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Taking into account the equality: Di{0,i) = gi{0,i),l < i < s, we derive the formula of a 
solution of the Cauchy problem (|T^.(|T6l) for initial data ( |T8l) 

s 

g,{t, Y) = 2t,(t, Y) H g,{0, z), s > 1, (19) 

i=l 

where ^sit,Y) is the s-order cumulantlf] of the groups of operators ([7]) of the von Neumann 
equations [28|[29] 

'^s{t,Y)= Yl (-1)'""'(|P|-1)! n (20) 

P:Y=\J^X, X,GP 

and we use accepted above notations. It should be emphasized that the structure of expansions 
fl20|) of the evolution operators 2ls(t, Y) means that they are determed by the cluster expansions 
of the groups of operators ([7]) of the von Neumann equation ([3])-(ll]). 

Thus, the cumulant nature of correlation operators induces the cumulant structure of a 
one-parametric mapping generated by solution (fT9l) and as we shall see in section 4, the same 
statement is true for the cumulant structure of the solution expansion of the BBGKY hierarchy. 

From (fT9|) it is clear that in case of absence of correlations in a system at initial instant the 
correlations generated by the dynamics of a system are completely governed by the cumulants 
(|2U]) of evolution operators ([7]). 

Hereafter we use the following notations: Yp = {{Xi}, . . . , {X\p\}) is a set, elements of 

which are |P| mutually disjoint subsets Xj C F = (1, . . . , s) of the partition F : Y = U^^^Xi, 
i.e. |lp| = |P|. In view of these notations we state that {Y} is the set consisting of one element 
F = (1, . . . , s) of the partition P (|P| = 1) and \{Y}\ = 1. 

We define the declasterization mapping 9 : Yp ^ Y, hj the following formula 

e{Yp) = Y. (21) 

For example, let X = (1, . . . , s + n), then for the set ({F}, X\Y) it holds 

e{{Y},X\Y) = X. 

Now we construct a solution of hierarchy (fT4|) for arbitrary initial data. Using solution ([6]) 
of the von Neumann equation ([3]) we rewrite formula (fT3!) as follows 

g,{t,Y) = Yl (-1)"''"'(|P| - ^y-^t n G\x.i{-t,X,)D\x.\{0,X^)- 

In view of expansion (Q at initial instant we have 

9s{t,Y)= J2 (-1)"''"'(|P|-1)!'5^ n ^l^^lH'^^) E U 9\zjO,Z,^). 

P-Y = U^X, X^CP P'^^«=U,, Z,^ Z,^ CP' 



■^Some properties of cumulants are considered in section 3.1. 
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Changing the summation indexes P and P' in this expression gives us 

9sit,Y)= E (-l)l^'l-^(|Pl-l)!x 

P : y = U,X, P':({Xi},...,{X|p|})=Ufe 

X n 0\eiz,)\i-t,eiz,))sf n m\io,x,). 

Finally we represent the solution of Cauchy problem 0141] . (1161) for arbitrary initial data in the 
following form 

gsit,Y)= Yl '^\p\it,{Xl},■■■,{XlPl})Sfl[g^xMX^)^ ^ > 1- (22) 

P:Y=[J^X, X,CP 

Here 2l|p|(t) is the |P|-order cumulant of groups of operators ([7]) defined by the formula 

2l|p,(t,{Xi},...,{X,p,})= (23) 
E (-I)I^'I^^(IP'I-I)! n S\eiz,)\i-t,0{Z,)), 

P':({Xi},...,{X|p|})=Ufc^fe ZkCP' 

where ^p'.({Xi} {X|p|})=Uj. possible partitions P' of the set ({Xi}, . . . , 

{X|p|}) into |P'| nonempty mutually disjoint subsets C ({Xi}, . . . , {X|p|}). 

We remark that for |P| > 2, the |P|-order cumulants 2t|p|(t) of groups of operators ([7]) of the 
von Neumann equations have similar structure in contrast to the first-order cumulant. Indeed, 
for |P| = 1 we have 

^^l{t,{l,...,s}) = gn{-t,l,...,s), 

and, for instance, in case of |P| = s, according to fl25]) . the s-order cumulant has structure (I2U|) . 
The simplest examples of correlation operators f l22|) are given by the expressions 

g,{t,l)=Qli{t,l)g,{0,l), 

g2{t, 1, 2) = 2ti(t, {1, 2})g2{0, 1, 2) + 2t2(t, 1, 2)StgiiO, l)(7i(0, 2), 

g^it, 1, 2, 3) = 2ti(t, {1, 2, 3})^3(0, 1,2,3) + 2l2(t, {2, 3}, l)Sig^{0, l)g2{0, 2, 3) + 

+2l2(t, {1,3}, 2)53±(7i(0, 2)(72(0, 1, 3) + 2l2(t, {1, 2}, 3)Sig^iO, 3)(72(0, 1, 2) + 

+2t3(t, 1, 2, 3)53±(7i(0, l)(7i(0, 2)(7i(0, 3). 

The validity of solution expansion f l2^ can be verified by straightforward differentiation by 
time variable (see section 3) and also in the following way. Taking into account the fact that 
the von Neumann hierarchy (fT^ is the evolution recurrence equations set, we can construct a 
solution of initial- value problem (fT^ .( !T6|) by integrating each equation of the hierarchy as the 
inhomogeneous von Neumann equation. For example, as a result of integrating the first two 
equations of hierarchy (ITl]) and using formula (^^, we obtain the following equalities 

g^{t,l) = g,i-t,l)g,{0,l), 

g2{t,l,2) = g2{-t,l,2)g2{0,l,2) + 

t 

J dhg2{-t + ti, 1, 2) ( - Af^;^{l, 2))Sig,{-tu l)gi{-h,2)g,{0, l)g,{0, 2). 



+ 
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In consequence of the fact that the operators commute with groups Gi{—ti, 1), ^i, 2), 
the second equahty takes the form 

g2{t,l,2) = g2i-t,l,2)g2{0,l,2) + 

t 

+ j dti6;2(-t + ti,l,2)(-Ari?(l,2))^i(-ti,l)6;i(-ti,2)52±(7i(0,l)(7i(0,2). 



Then for the second term on the right-hand side of this equation an analog of the Duhamel 
equation holds 

t 

j rfti6;2(-t + ti,l,2)(-A/;i?(l,2))6;i(-ti,l)6?i(-ti,2) = (24) 



t 

/d 
dti—{g2{ti, 1, 2)g,{-h, i)g,{~ti,2)) = 



= g^i-t, 1, 2) - gii-t, i)gi{-t, 2) = 2i2(t, 1, 2), 

where 2l2(t) is the second-order cumulant of groups of operators ([7j) defined by formula (120|) . 
For s > 2 the solution of the Cauchy problem (fT4l) .( fT6l) . which is constructed by iterations, is 
represented by expansions ( |22l) as a consequence of transformations similar to an analog of the 
Duhamel equation ( HM . 

We note, that in case of initial data f lTSj) solution f|T9|) of the Cauchy problem f|T^ . f|T6|) of 
the von Neumann hierarchy may be rewritten in another representation. For n = 1, we have 

gi{t,l) = Qli{t,l)9i{0,l). 

Then, within the context of the definition of the first-order cumulant, 2li(t), and the dual 
group of operators 2li(— t), we express the correlation operators gs(i), s > 2, in terms of the 
one-particle correlation operator gi{t) using formula ( !T9l) . Hence for s > 2 formula (fTOll is 
represented in the form 



gs{t,Y) = Qis{t,Y)S^ Yl gi{t,i), s > 2, 

1=1 

where Qls{t,Y) is the s-order cumulant (120|) of the scattering operators 

s 

gt{Y) = gs{-t,Y)l[gi{t,k), s>i. 



k=l 



The generator of the scattering operator gt(Y) is determined by the operator 

pt{Y)U=o = -i2 E Ar(?(^i,...,^,), 

k=2 ji<...<ifc = l 

where the operator Af^^^^ acts in £q(?/^) C £,^{1-Lf) according to formula (|T5l) . 
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2.5 Correlation operators of particle clusters 

Further we consider a more general notion, namely, the correlation operators of particle clusters 
that describe the correlations between clusters of particles. For a system consisting of particles 
and of one particle cluster {Y} = {1, . . . ,s} we introduce the sequence of operators g^^\t) = 
{gi+oit, {Y}), . . .,gi+n{t, {Y},X\ Y), ...) e £H©^=o'^^+n) that are defined by the following 
generalized cluster expansions of the density operators D{t) = (1, 1), . . . , Ds+nii) 1, . . . , s+ 



n),...)G£H-^^ 



HI 

Ds+n(t, X) = g^+r,{t, {Y},X\Y)+ (25) 

+ E Sf^^\[g\xS^X{), n>0, 

P:({y},x\y) = Ui-Y., x,cp 

|p| > 1 

where X = (1, . . . , s + n), X\Y = (s + 1, . . . , s + n) and Z^P: ({y},x\y)=u ^i, |p|>i ^^"^ 
over all possible partitions P of the set \ Y) into |P| > 1 nonempty mutually disjoint 

subsets Xi C {{Y},X\Y). 

Relations between correlation operators of particle clusters g^^\t) = (gi+o(t, {Y}), . . . , gi+n{t, 
{Y},X \Y),. . .) G £^(©J^o^s+n) correlation operators of particles (fT3|l is given by the 
equality 

g,+r,it,{Y},X\Y)= (26) 

E (-i)'"-'(ipi-i)!'5^+. n E n 

In particular case n = 0, i.e. the correlation operator of a cluster of |y | particles, these relations 
take the form 

g.^oitAY}) = St E Uaix.lii^X,). 

P:Y=[J,X, X,CP 

From relations fl26|) using hierarchy f|T^ we derive the von Neumann hierarchy for the cor- 
relation operators of particle clusters 

"^-gi+nit, {Y}, X\Y) = -A/;+„(X)(7i+„(t, {Y},X\Y)+ (27) 



dV 



(E \z^\) 



P : {{Y}, X\Y) = \J^X,, ZiCe(Xi), Z|p| C e(X|p| ), 

|P| > 1 Zi^9 Z|p| 5^ 

Z\P\))St,n n 

The formal solution of the von Neumann hierarchy for correlation operators of clusters of 
particles is given by the following expansion 

^i+„(t,{F},X\F)= (28) 

J2 "^IPlit, {0{Xi)},...,{0{X\P\)})Sf_,^ II g\xA^,X,), n>0, 
P:{{y},x\y)=ux, x,cp 
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where 2l|p|(t) is the |P|-order cumulant defined by ( l23l) . 

For the sequence of correlation operators of clusters of particles g^'^\t) = {gi^o{t, {Y}), gi+i{t, 
{Y}, s + 1), . . . , gi+ni{y}, s + 1, . . . , s + n), . . .) the initial data satisfying a chaos property, i.e. 
in case of absence of correlations between particles and between particle cluster and particles, 
is the sequence 

g^'\0) = ig,^o{0,{Y}),0,...). (29) 

The corresponding solution of the initial- value problem of the von Neumann hierarchy is given 
by the expansion 

<7i+„(t,{y},x\y) = 2ti+„(t,{F},x\y)5,=^+„ J] ^1(0,2), (30) 

ie{{Y},X\Y) 

where 2li+n(t) is the (n + l)-order cumulant defined by formula (!23|) and has the form 

Qi,^^{t,{Y},X\Y) = Yl i-^f^-'m - 1)! n G\oixd-t^e{X,)). (31) 

P:({y},x\y)=U;X, x,cp 

Thus, correlations created during the evolution of a system are described by formula fl30l) and 
defined by the corresponding order cumulant fl3Tl) . 



3 The Cauchy problem of the von Neumann hierarchy 
in the spaces of sequences of trace-class operators 

3.1 Some properties of a solution 

The intrinsic properties of constructed solution (l22l) are generated by the properties of cumu- 
lants fl251) of groups of operators of the von Neumann equations. 

At initial time t = solution (!22|) satisfies initial condition (fT6|) . Indeed, according to 
definition ([7j), i.e. that ^n(O) = / is a unit operator, and in view of the identity 

n 

E (-i)'""'(ipi - 1)! = Y.(-ir-Mn, m - ly. = 6^^„ (32) 

P:{l,...,n)=U,X, k=l 

where s(n, k) are the Stirling numbers of the second kind and ^^^i is a Kronecker symbol, for 
|P| > 2 we have 

2l|p|(0, {X,}, . . . , {X|p|}) = (-1)"'''"'(|P'I - 1)!^ = 0. 

P':({Xl},...,{X|p|})=UfeZfe 

Let us consider some properties of cumulants fl20|) of groups of operators ([7]) introduced in 
previous section. 

Since for a system of non-interacting particles it holds 

n 

Qni-t, 1, . . . , n) = JJ Gii-t, z), 

i=l 
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then taking into consideration identity ( l32l) . for n > 2 we obtain 

2i„(f,i,...,n)= (-i)""'-'(iPi-i)! n n^iH'^^)= 

P:(l,...,n)=U,X, X,<ZPh=l 
n n 

= Y^(-lf-^sin,k)ik-l)\l[g,{-t,t) = 0. 

k=l 1=1 

The infinitesimal generator of the first-order cumulant is defined by the following limit in 
the sense of the norm convergence in the spaces £^('H^) 

lim i(2ti(t, {!,..., n}) - /)/„, = (-Ar„/„)(1, ...,n), 

where the operator (— A/'„) is defined by formula for /„ G £o(^n) '^^i.'^t)- 

In general case, i.e. n > 2, the generator of the ri-order cumulant is given by the operator 

lim 1, . . . , n)U = (-Ari!:Vn)(l, ...,n), (33) 

defined by formula ffT^ on £,1{'H^) C SMTi^) and the limit exists in the sense of the norm 
convergence in the spaces £^('H^). To prove the validity of statement (133|1 we use definition 
(!20|) and equality ([8]) and represent the limit in left hand side of (133!) in the explicit form 

limi2l„(t)/„ = limi (-I)I^I-'(IPI-I)! \{ <3\z,\{-h Z,) U = (34) 

^ P:(l,...,n)=Ufc2fe ZfcCP 

P:(l,...,n,)=Ufc2fc ^fcCP 

n 

= E (-i)'"^'(ip|-i)'E(--^i(^))/" + 

P:(l,...,n)=Ufc^fc i=l 

+ E (-i)'"-^(ip|-i)!E E E i-^i:lin^---^3n))fn, 

P:(l,...,n)=UfeZfe ^feCP r=2 ji<...<j,=r 

where Afi{i)fn = \{K{i)fn — fnK{i)). Taking into account identity ( 132|) and the fact that 
n > 2, for the first term of formula fl34l) it holds 



E (-i)'"-^(ip|-i)!E(--^i«) = o- 

P:(l,...,n)=Ufe2fc i=l 

Applying identity ( !32l) to the result of rearrangements in the second term of ( l34l) . we finally 
obtain 

E (-iri-^(ip|-i)!EE E (-<^(ji,---,^.)) = 

P: (l,...,")=Ufe Zk ZkCP r=2 ji<...<>=r- 
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n~\Z\+l 

= E(--^i2?^(^)) E (-i)"i-^(|p|-i)!s(n-izi + i,fc) = 

ZCY k=l 

= E(--^2f'^(^)) 5n,izi = -^f!^il, . . . 

ZcY 

We emphasize that the operator (— A/'^"^) is defined by a n-body interaction potential ([2]). 

Similarly to the proof of statement fl33l) for the |P|-order, |P| > 2, cumulant fl23l) of particle 
clusters for /„ G iigl^n) -^H^n) "^^ obtain in the sense of the norm convergence in the 
spaces £,^(T-L^) 

liinl2l|p|(t,{Xi},...,{X|p|})/„= (35) 

t-s-O t 

E (-i)"'''"'(ip'i - 1)! E i-^oiz.mz,)))/^ = 

P':({Xi},...,{X|p|})=U,^fc ZfeCP' 

|P| 

= E E (--^in^ (^l,...,^|P|))/„. 

-^1 C Xi, ^\P\ ^|P|' 

^1 7^ Z|P| 7^ 

3.2 The group of nonlinear operators generated by the von Neu- 
mann hierarchy 

In the spaces constructed solution (1221) generates a group of nonlinear operators. The 

properties of this group are described by the following theorem. 

Theorem 1. For f e £H-^-h) 

mapping 

i ^ = E 2t|p,(t, {XJ, . . . , {X|p|})5„± n MX.) (36) 

P:(l,...,n)=U,X, X,CP 

is a group of nonlinear operators of class Cq. In the subspaces £,1(71^) C 2^(1-1^), ny 1, the 
infinitesimal generator 9T(-) of group is defined by the operator 

(0I(/))Jl,...,n) = -A4/n+ (37) 
+ E E ••• E {-MF'^\z„...,Z^^^))S^l[f\x.\{X,), 

P:(l,...,n) = U,^,:, ^iCXi, Z|p|CX|p|, ^iCP 

|P| 7^ 1 ^1 7^ Z|P| 

Proof. Mapping fl36l) is defined for /„ G -2^(7/^), -n, > 1, and the following inequality holds 

where c = max(c, 1) and c = maxp. (i ..._„)=u.Xi ,i ± • Indeed, since for /„ G 2^(71^) it 

holds (Proposition 1) 

Tri,...,„|^„(-t)/„| = ||/«||£i(^±), n>l, 
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hence taking into account the inequahty 

W^n n •^l^^lll 1 ± - n ll-^l^'lll , ± ' 

we obtain 

P:(l,...,n)=UiXi p':({Xi},...,{X|p|})=UfeZfc ^>CP 

< 5: clPIX:s(|P|,^)(fc-l)!< E clP'E^'''-^ = 

P:(l,...,n)=U^> fc=l P:(l,...,n)=U,X> k=l 

n m 

m=l fc=l 

where s(|P|,A;) are the Stirhng numbers of the second kind. Thus, (2l((/))^ G ii^('H^) for 
arbitrary t G and n > 1. 

We now prove the group property of a one-parametric family of the nonlinear operators 
that are defined by fl36|) . i.e. 

2tt.(2ti,(/)) =2l*,(2t,,(/)) =2li,+i,(/). 
For fn G £^(7/^), n > 1, and for arbitrary ti, t2 E , according to definition (!23|) . it holds 
(21*. {%,if)))Jl, . . . , n) = ^ 2t|p|(ti, {Xi}, . . . , {X|p|}) X 

P:(l,...,n)=UX, 

X n E %.i(^2,{^ij,...,{v.ij)'5^ n /i^^i(^'.)= 

= E E (-i)"'''"^(ip'i-i)! n wHi'^w.))x 

P:(l,...,n)=UX. p':({Xi},...,{X|p|})=U,Q, Q,CP' 

xJI E E (-1)I^^I-^(|P:|-1)!x 

X,CP P,:X,=U^Z,, P:;({ZiJ,...,{Z,p^|J)=U,^R,, 

where {{Zi.}, . . . , {Z|p.|-}) is a set, elements of which are |Pi| subsets Zi- C Xi of the partition 
Pj : = Uj.Zi-. Collecting similar terms of products of operators /„, n > 1, and taking into 
account the group property of operators Qn{~t), n > 1, we establish 

(2l*.(2li,(/)))„(l,...,n) = 

= E E (-i)"''^^(|p'i-i)! n ^I^W^)lHi-^2,^(g.))x 

P:(l,...,n)=UX, P':({Xi},...,{X|p|})=U,Qfc QfcCP' 

n /i^.ira= E %i(ti+t2,{xi},...,{x,p,})5„± n /i^dra= 

X,CP P:(l,...,n)=U,Xi X,CP 

= (2l*.+*.(/))„(l,...,n). 
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The strong continuity property of the group 2tt(-) over the parameter t G is the straight- 
forward consequence of the strong continuity of groups of the von Neumann equations |23j . 
Indeed, according to identity (!32l) the following equality is valid 

E E (-l)"''''(|P'|-l)!'5„^ n /|x.|TO = /„(l,...,n). 

P:(l,...,n)=UX, P':{{Xi},...,{X|p|})=Ufc^fe ^«CP 

Therefore, for /„ G S^HU^) C S.\nt), n > 1, it holds 

P:(l,...,n)=UX. p':({Xi},...,{X|p|})=Ufe^fc ZfcCP' 
^'CP P:(l,...,n)=U,X,; P':({Xi},...,{X|p|})=Ufc^fe 

iim|| n Gmz,)\{-t,9{z,))s^ n n /i^.k^^Lh^^)- 



X 



ZfcCP' ^>CP X,CP 



In view of the fact that the mapping Qn{—t) is strong continuous group, i.e. in the sense of the 
norm convergence in £^(7/^) there exists the limit 

lim(6;„(-t)/„ - = 0, 
this implyies that for mutually disjoint subsets Xi C Y, the following equality is also valid 

ZfeCP' 

For fn e S^liU^) C 2}{U^), we finally establish 

li-||(2l*(/))„-/n|U^.) = 0. 

To construct the generator OX(-) of strong continuous group we first differentiate it in the 
sense of the pointwise convergence in the spaces i.e. for arbitrary ifjn G V{Hn) C T-L^. 

Taking into consideration the fact that for /„ G £0(7/^) C V{^{-)n) equalities ([8]) and (135!) 
hold, for group (!36|l we derive 

limi((2l,(/))„-/„)^^„= (39) 

= ^^0l^ 5Z %|(t,{^l},...,{X|p|})5„± n /|Xd™-/n)V^n = 

~^ P:(l,...,n)=U,X, X,CP 

= lini ^ (2ti(t, {1, . . . , - + 

+ E liml2t,p,(t,{Xi},...,{X|p|})5;^ J] /|x.i(X.)^n = 

P:(l,...,n) = U.X,, ^ ^>CP 
|P| > 1 

= {-Kf„)A+ E (-A/'^^H{Xj,...,{X|p|}))5„± n MXi)i^n, 

P : (l,...,n) = U,X,;, X,CP 
|P| > 1 
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where the operator A/""^* is given by the formula 

Ar'°*({Xi},...,{X|p|})^ J2 ■■■ E MF'^\z,,...,Z^p^). (40) 

Zi C Xi, Z|p| C ^|p|, 

Hence in view of equality for /„ G £o(^^) ^ ^^(^(On) C 2.\n^), n > 1, in the sense of 
the norm convergence in ii^('Hj), we finally establish 

lim||^((W))„-/„)-(^(/))J|,,(^,) = 0, 
where is given by formula ( l37j) . □ 

3.3 The existence and uniqueness theorem 

Let Y = (1, . . . , s) is above accepted notation. For abstract initial- value problem (ini) . (ITB]) in 
the spaces £}{J-'^) of sequences of trace class operators the following theorem is true. 

Theorem 2. The solution of initial-value problem [l4\), fT^) of the von Neumann hierarchy is 
determined by the following expansion 

gs{t,Y)= E '^\p\{t,{Xi},...,{X\P\})St H gix^i{0,X,), s>l, 

P:Y=[J^X, X,CP 

where 2t|p|(t) is the \P\-order cumulant of the groups of operators ^ defined by formula ^2&^) . 

For QnifS) G i2o('^n) -^^(^n) strong (classical) solution and for arbitrary initial 

data 5'„(0) G S^il-L^) one is a weak (generalized) solution. 

Proof. The Theorem 1 implyies that for initial data 5'„(0) G £.\{1-V^) C 2}{1-V^), n>l, sequence 
f l22|) is a strong solution of initial- value problem f|T^ .f lT6|) . 

Let us prove that in case of arbitrary initial data 5^(0) G S}{J^^) expansion fl22|) is a weak 
solution of the initial- value problem of the von Neumann hierarchy ( !T^ . With this purpose we 
introduce the functional 

{fs^gs{t)) = Tn,...,, . . . , s) g,{t, 1, . . . , s), (41) 

where fs G £o(^^) is degenerate bounded operator with infinitely times differentiable kernel 
with compact support and the operator gsif) is defined by expansion ( 122|) . According to estimate 
dSH]), for /, G S.Q{Hf), functional (jH]) exists. 

Using expansion fl2^ . we transform functional f HT]) as follows 

{fs,gs{t))=Tn_s E /.%l(i,{^i},---,{^|P|})5f n ^?|x.|(0,X,)= (42) 
= Trw E E (-l)l^'l-^(|P'|-l)!x 

P: y=U, X. P': ({Xi},...,{X|p|})=Ufe 

X n ^ie(z,)i(t,^(^fc))/.5± n 9\x.\i^^Xi)^ 

ZfeCP' ^»CP 
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where the group of operators G\e{Zk)\(i) is adjoint to the group Q\9(Zk)\{~t) iii the sense of 
functional ( HTl) . 

For gn{0) G ii^('H^) and /„ G ilo('H^) within the context of the Proposition 1 [22] the 
following equality holds 

lim Tri,...,„ {^AGn{t + At)fn - 6^n(t)/„) - Gnit)N'nfn)gn{0) = 0, 
At->-0 Ac 

and, therefore we have 

ZfeCP' ZfeCP' 



X,CP 

= Tri„„,, 5^ A/],(z,)|(^(^0)/.(n n G\eiz,)i{-t,9iZ,))St H g^x^^iO, X,). 

ZiCP' ZkCP' -^iCP 

For fs G iio(^^) and bounded interaction potentials the limit functional exists. Then, taking 
into account that symmetrization and anti-symmetrization operators ([1]) are integrals of motion, 
the derivative over time variable of functional (HTj) represented by (142!) is given by 



IP 1-1 X 



I (A, 9s{t)) = [Afsfs, gs{t)) + Tt,_s J2 E (-1) 

P: Y = [j^X,, P':({Xi},...,{X|p|})=Ufe2fe 

|P| > 1 

x(|p'|-i)! E ^\eizomzi))fs{Y) n gm,)i{-t,e{z,))st n ^?|x.|(0,^.) = 

Z,CP' Z,CP' ^»CP 

= (AC/s, +Tri,.„,, Yl Ar-*({Xi},...,{X|p|})/,(y)5± n 

P:y = U^», X»CP 
|p|>i 

where the operator Af™^{{Xi}, . . . , {X\p\}) is defined by formulae fHOj) and f|T5|) . 

Thus, the result of the differentiation over time variable of functional pT]) takes the form 



j^Us.dsit)) = {Kfs,gs{t))+ (43) 

|p|>i 

Equality (H3!) means that for arbitrary initial data (7,1(0) G £^(7^^), n > 1, a weak solution 
of the initial- value problem of the von Neumann hierarchy (fT4|) . (fT6|l is determined by formula 
(121. □ 



The statement of this theorem also holds for the correlation operators of clusters of particles 
introduced in section 2.5. 
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4 On applications of the von Neumann hierarchy 

4.1 Functionals of the mean- value and the dispersion of observables 

Observables of the Bose and Fermi many-particle systems are sequences A = {Aq, Ai, . . . , 
An, . . .) of self-adjoint operators defined on the Fock spaces J-"^ and Aq G C. We consider 
observables as elements of the spaces of sequences of bounded operators S.{T^) with an operator 
norm [^UlES]- The average values of observables (mean values of observables) are defined by 
the positive continuous linear functional on the spaces 2.{J^^) 

(A) = {I, D)-' V - Tri,...,„ An Dn, (44) 

n=0 

where Tri ... „ are the partial traces over 1, . . . , n particles, (/, D) = Xl^o'^rTri,. ..,n-D„ is a nor- 
malizing factor (grand canonical partition function) and D = [1, Di, . . . , Dn, . . .) is a sequence 
of self-adjoint positive density operators defined on the Fock spaces J-'^ that describes the 
states of a quantum system of non-fixed number of particles. For D G il^(J-'^) and A G 
average value functional fHijl exists. It determines a duality between observables and states. 

We remark that in case of a system of fixed number of particles the observables and 
states are one-component sequences, respectively, A*^^) = (0, . . . , 0, A^v, 0, . . .), D^^^ = (0, . . . , 
0, Dtv, 0, . . .), therefore, formula (jH]) for an average value reduces to 

(AW) = (Tri,...,jv/^jv)"'Tri,...,^A^D^. 

In the framework of the description of states by correlation operators f|T3|) the average values, 
for example, of additive-type observables A^^^ = (0, ai, . . . , ^"^j^ai(z), . . .) are given by the 
functional 

oo ^ 

(A«)(t) = y^-, Tri,...,i+„ ai(l)f/i+„(t, 1, 2, . . . , 1 + n), (45) 

n=0 

or in general case of s-particle observables A^^^ = (0, . . . , 0, 0^(1, . . . , s), . . . , Ylh< <is=i ^sihi ■ ■ ■ , 
is), ■ ■ .) by the functional 

1 °° 1 

= - V - Tri,...,,+„ a,(l, . . . , s)gi+n{t, {1, . . . ,s},s + 1, . . . ,s + n), (46) 
s! ^ — ^ nl 

n=0 

where the correlation operator gi-^-n{t) is defined by expansion (EH]). For A(^) G £(7"^) and 
g^'^ G functional 1^ exists. 

In view of the definition of average value functional f l46p . the dispersion of an additive-type 
observable is defined by the functional 

oo ^ 

((^« - (^«)(t))^)(t) = J2^_ Tri,...,i+„ (a?(l) - (A«)^(t))^i+.(t, 1, . . . , 1 + n) + 

n=0 

oo ^ 

+ ^ Tri,...,2+„ ai(l)ai(2)(72+„(t, 1, . . . , 2 + n), (47) 
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where {A^^^)(t) is determined by expression gHD. For ^(1) G S.{J^n) and g E ^^{J^u) functional 
(HZD exists. 

To justify definition (146|) . i.e. to derive functional (146|) from ( H4|) . we introduce necessary 
notions and formulate some equalities. 

For arbitrary / = (/o, /i, ...,/„,.. .) G and Y = (1, . . . , s) we define the linear 

mapping : / — )■ ^y/ by the formula 

(f y/)„ = fwi+niY, |y I + 1, . . . , |y| + n), n > 0. (48) 
For the set {Y} defined above in section 2.4 we have respectively 

i^{Y}f)n = fi+ni{Y}, s+l,...,s + n), Tl > 0. (49) 
On sequences dyf and Oy/ we introduce the ©-product 

(0y/®3y,7)|x|(X) = 5^5±^u^,u^)|/|^| + |y|(r,Z)7[x\Z| + |y'|(r,X\Z), (50) 

where Sp^ are operators defined by formula ([I]), 6' is the declasterization mapping defined by 
fl2T|) . X, y, y are the sets, elements of which characterize clusters of particles, J2zcx 
sum over all subsets Z of the set X. In particular case {Y = 0, F' = 0) definition f l50|) reduces 
to ([10]). 

For / = (0, /i, ...,/„,.. .), /„ G £^(7/^), according to definitions of mappings ( ITTj) and (jlSD 
the following equality holds 

0{y}Exp^/ = Exp^/ ® d{Y}f, (51) 

and for mapping (148|) correspondingly, 

c)yExp^/ = Exp^/ ® J] OxJ®...®fX|p|/, (52) 

where ^p.y=u,x, ^^m over all possible partitions P of the set Y = (1, . . . ,s) into |P| 

nonempty mutually disjoint subsets Xi C Y. 

For /„ G £^(?^^) an analogue of the annihilation operator is defined by the formula 

(a/),(l, ...,s) = Tr,+i/,+i(l, ...,s,s+l), (53) 

and, therefore 

oo ^ 

(eV)s(l, • • • , s) = -rTrs+i,...,s+n/s+n(l, • • • , s + n), s > 0. 

n=0 

According to definitions fl53|) and fl50l) . for sequences /, / G I1^(J-'^) the equality holds 

(e"(/®7))o = (eV)o(eV)o. (54) 
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In terms of mappings ( HHl) and ( H9|l generalized cluster expansions ( 125|) take the form 

OyD(t) = r){y}Exp^ g{t), (55) 
We now deduce functional fj46l) from the definition of functional (jH]), i.e. 

= i(e"D(0))o^Trya,(y)(e"?)yD(t))o, 

where Try = Tri_ Indeed, using generalized cluster expansions (l55ll and as a consequence 
of equalities ( 15T|) and ( 15^ . we find 

(e°D(0))o^(e^OyD(t))o = (e"D(0))o ^(e^ 0{y}Exp^ ^7(t))o = 

= (e"D(0))oi(e"Exp^(7(t) ®0{y|(?(t))o = (e^D(0))o i(e'^Exp^^(t))o(e''?)|y|^(t))o. 

Observing that Exp(g(7(t) = D{t) and (e°D(0))o = (e"L)(t))o, as a final result we derive repre- 
sentation (H6l) . 

We note that the possibility of the description of states of systems in the framework of 
correlations arises naturally as a result of dividing of the series in expression (jUj) by the 
normalizing factor series. Indeed, the solution of generalized cluster expansions fl55|) has the 
form 

^{Y}9{t) = 0{y}^9Lng,D(t). 
Here for f\Yp\{Yp) e S^CH^,) the mapping is introduced by 

W)|y,|(yp) = /|,(y,)|(e(rp)) = /|y|(y), 

where 9 is the declasterization mapping (12T]) defined on sets Ip = ({Xi}, . . . , {X|p|}) consisting 
of clusters of particles. Then in terms of a sequence of the density operators D{t) the following 
representation of functional (jUj) holds 

= iTrya,(y)(e"0{y}^9Lng,D(f))o. 

It should be emphasized that correlation operators that belong to the spaces £}{J-'^) describe 
only finitely many particles, i.e. systems with finite average number of particles. Indeed, for 
the observable of number of particles N = (0, /, 21, . . . , nl, . . .) functional fHSj) has the form 

oo ^ 

(N) {t) = y^-r TTl,...,l+n 9l+n{t, 1, . . . , n + 1), 

n=0 

and, for example, in case of initial data satisfying chaos property (fT9|) . if || (71(0) ||£i(^) < e~^, 
we get 

00 

\{N){t)\ < e||^i(0)||£i(^) 

n=0 
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4.2 Marginal density operators 

In view of functional (146|) we define the marginal density operators Fs{t) in terms of correlation 
operators of clusters of particles g^'^\t) = (5'i+o(^) {^})) ■ ■ ■ i Oi+niii {Y}, s + 1, . . . , s + n), . . .) 
by the expansion 

oo ^ 

s+l,...,s+n 9l+n (t,{y},s + l,...,s + n), s>l, (56) 



n 

ra=0 



where 5fi+„(t, {Y}^ s + 1, . . . , s + n), n > 0, is solution ( l28l) of the von Neumann hierarchy (127|) . 
Then the average value functional of s-particle observable f H6l) is given by the formula [l£ 



= 1 Tn,...,, a,(l, . . . , s)F,(t, 1, . . . , s). 
s\ 

The initial- value problem for the marginal (s-particle) density operators (!56|) has the follow- 
ing form (BBGKY hierarchy of the Bose and Fermi many-particle systems) [30] 

j^F,{t,Y) = -Ns{Y)F,{t,Y) + 

oo ^ 

+ E ;^Tr,+i,...,,+„ (-A/;2f'^"^)(^,^ + l,---,^ + ^)^s+n(t), (57) 

n=l ■ Z CY, 

Fs{t)\t=o=Fs{Q), s>l. (58) 

We remark, that as follows from f lH|) in the framework of the description of states by the 
density operators D = {l,Di{t), . . . , Dn{t), . . .) the marginal density operators Fs{t), s > 1, are 
defined by the well-known formula (the case of the non-equilibrium grand canonical ensemble 

oo ^ 

Fs{t, l,...,s) = {I, D{0))-' V -TTs+i,...,s+nDs+n{t, l,...,s + n), 

or in terms of definitions ( HHj) and fl53|) . we rewrite this formula in the following form 

i^y|(t,F) = (e"D(0))oi(e"OyD(t))o. 

Then as above owing to generalized cluster expansions fl55|) and the validity of equalities f lSTj) 
and flM|) . we deduce expansion fl56l) 

F|y|(t,y) = (e"c){y}^7(t))0. 

Thus, along with the definition in the framework of the non-equilibrium grand canonical ensem- 
ble the marginal density operators can be defined in the framework of dynamics of correlations 
that allows to give the rigorous meaning of the states for more general classes of operators than 
trace class operators. 

In section 4.4 on basis of the solution of the von Neumann hierarchy (!27ll the solution of 
Cauchy problem ( l57j) - (!58l) of the BBGKY hierarchy for the Bose and Fermi many-particle 
systems is constructed. 
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4.3 Marginal correlation operators 

In view of ( 147|) we introduce the marginal correlation operators by the series 

oo ^ 

Gs{t,l, . . . ,s) = y^ —Tis+i,...,s+n gs+n{tA, ■ ■ ■ ,s + n), s > 1, (59) 

n=0 

where gs+n(t, 1, . . . , s + n), n > 0, is solution ( !22l) of the von Neumann hierarchy ( fT4l) . Thus, 
dispersion functional ( 1471) represents as follows 

- = Tn (a2(l) - (A«)2(t))G'i(t,l) + Tri,2ai(l)ai(2) G2(t,l,2). 

We note, that the evolution of marginal correlation operators ( 159|) is governed by the initial- 
value problem of the nonlinear BBGKY hierarchy [19] , [26] . 

Marginal correlation operators (l59l) are interpreted as the operators that describe correlations 
of the marginal states (156|1 . In fact, it holds 

Gs{t,Y)= (-1)"''"'(|P|-1)''5„^ n ^^1' (60) 

P:Y = \J,X, X,CP 

i.e. the marginal correlation operators are the cumulants (semi-invariants) of the marginal 
density operators. In terms of mapping ( JT2l) sequence (|60l) can be rewritten in the form 



G{t)=hn^F{t), 

where F(t) = 1 + (0, Fi{t), . . . , Fn{t), . . .). As consequence in terms of the marginal density 
operators dispersion functional ( H71) takes the known form [T9| 



((A« - {A'^'^) f){t) = Tn (alii) - (A«)2(t))Fi(t, 1) + 
+ Tri,2 ai(l)ai(2) {F^it, 1, 2) - 52±Fi(t, l)Fi(t, 2)) . 

Let us prove the validity of relations fl60|) . Consider expression (OyExpg,G'(t))o. In view of 
expansion (15^ . using definition flTT]) of the mapping Exp^, and equality f lM|) . we obtain 

(c)yExpgG(t))o = (OyExp^e"(7(t))o = 

P.Y=\J,X, X,CP P:Y=\J,X, 

Then according to equahties (!54l) and (!52|) . we have 

(e"Exp^^(t))o '(e"Exp^(7(t))o(e° Ox,^7(t) ® • • • ® ^x,^,g{t))o = 

P:Y=\J^X, 

= (e"Exp^(?(t))o-i(e^(Exp^(?(t) ® J] ^x,^7(t) ® • . . Ox|p|^(t)))o = 

P: 'i^=U. X, 

= (e"Exp^(7(t))oi(e»r)yExp^(/(t))o. 
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Observing that from relations ( l26l) it follows that 

and using equalities f lSTj) and f l5^ . we derive 

(e"Exp^^(t))o^(e"OyExpg,^(t))o = (e"Exp^(7(t))o i(e"c){y|Exp^(7(t))o = 
= (e"Exp^^7(t))o"'(e^Exp^(7(f) ® 0|y}(7(t))o = {e'^{Y}g{t))o. 

Therefore, taking into account definition we finally establish the validity of the cluster 
expansions 

F(t) = Exp,^G'(t). 

The solutions of obtained recursive relations in component- wise form are expansions fl60l) . 

4.4 On solutions of the BBGKY hierarchies in case of chaos initial 
data 

Let us apply introduced expansions of marginal correlation operators (159!) and marginal density 
operators (156!) to construct solutions of the Cauchy problems of the corresponding BBGKY 
hierarchies for one physically motivated example of initial data, namely initial data satisfying 
chaos property (ITS]) . 

In terms of marginal correlation operators fl59l) a chaos property means that 

a(0,l,...,s) = (7i(0,l)5,,i, s>l, (61) 

where 5s^i is a Kronecker symbol. Taking into account the structure of solution fll9p of the initial- 
value problem of the von Neumann hierarchy f|T^ in case of initial data f|T8l) . for expansion 
fl59|) we obtain 

oo ^ s+n 
Gsit,l,...,s) = ^ — Tls+l,...,s+n^s+nitA,---,S + n)S^_^^ Jj5'i(0,i), 

n=0 i=l 

where 2l<j+„(t) is the (s + n)-order cumulant (1201) . In consequence of (I6T1) we finally derive 

oo ^ s+n 

Gs{t, 1, • • • , = 5^ ^ Tr,+i,...,,+„ 2t,+„(t, 1, . . . , s + n)5,=^+„ J] ^(0, z), s > 1. (62) 

n=0 ■ i=l 

Since the estimate holds 

P„(t)/l|£i(^±) < n!e"||/||^i(^±), 

series (!62|) converges, if ||G'i(0)||£i(-^) < e^^. 

Thus, the cumulant structure of solution ( 122|) of the von Neumann hierarchy (fT4|) induces 
the cumulant structure of solution expansion (162 p of the initial-value problem of the nonlinear 
BBGKY hierarchy for marginal correlation operators. 
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In terms of marginal density operators ( l56l) chaos property ( l29l) means that 

F,(0,F) = (7i(0,{y}), s>l, (63) 

where Y = (1, . . . , s) and {F} is a notation introduced in section 2.4. According to formula f l30|) 
for a solution of the initial-value problem of the von Neumann hierarchy ( [T4l) from definition 
( 156|) we obtain 

F,(t, = 5^ Tr,+i,...,,+„ 2li+„(t, {y}, X \ Y) 5± „ J] ^^(0' 

n=0 ' i&{{Y},X\Y) 

where 2li+„(t) is (n + l)-order cumulant (l3Ti) . From definition (1561) in view of equality (1631) . i.e. 
fl'i(0,«) = Ilj(^0(i)Fii^J)^ we finally derive |30j 

oo ^ s+n 

F,(t,F) = ^-Tr,+i,„„,+„2ti+„(t,{y},X\r)5±„ H^i^^'^)' ^ ^ (^4) 

n=0 ■ i=l 

We point out that in case of chaos initial data obtained solution expansion (!64l) differs from 
solution expansion ( l62l) of marginal correlation operators only by the order of the cumulants of 
the groups of operators of the von Neumann equations. 

For arbitrary initial data F(0) G S.l{J^^) = ®'^=o'^"' i'^t) ^ if a > e, for t G M there exists 
a unique solution of the initial value problem f l57|) -f l58|) given by the series 

oo ^ 

F,(t,y) = V- Tr,+i,...,,+„,2li+„(t,{r},X\y)F,+„(0,X), s>l, (65) 

n=0 

where 2ti+„(t) is {n + l)-order cumulant (j3Tl) of groups of operators ([7]). For initial data 
F(0) G -Cq oI-^-h) ^ -^al-^-w) it is a strong solution and for arbitrary initial data from the spaces 
^lyi-^u) it is ^ weak solution. 

It should be emphasized that the solution expansions fl65l) for the Fermi and Bose many- 
particle systems has the same structure, as in case of the Maxwell-Boltzmann statistics [28] 
which caused by the fact that symmetrization and anti-symmetrization operators ([1]) are inte- 
grals of motion. 

We remark that a solution of initial- value problem (IFr|) -f l55|) is usually represented by the 
perturbation (iteration) series [2], [3], [7], [8], [13], [32]. For certain classes of interaction po- 
tentials on suitable subspaces of £^(J-'^) for group ([7]) and cumulants (!3T!) the analogues of the 
Duhamel formula hold [33| and as a result expansion (165|1 reduces to the iteration series of the 
BBGKY hierarchy. For classical systems of particles the first few terms of cumulant solution 
expansion for the BBGKY hierarchy were established in [26], j31j . 

5 Conclusions 

In the paper it was considered the origin of the microscopic description of non-equilibrium 
correlations of Bose and Fermi many-particle systems. For the correlation operators that give 
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an alternative description of the quantum state evolution of Bose and Fermi many-particle 
systems, the von Neumann hierarchy of nonlinear evolution equations ( fT^ was introduced. 

The properties of group of nonlinear operators generated by the von Neumann hierarchy 
were examined (Theorem 1). It gave us the possibility to construct solution (!22|) of the Cauchy 
problem of this hierarchy in the corresponding spaces of sequences of trace class operators 
(Theorem 2). 

In particular, it was established that in case of absence of correlations in the system at 
initial time the correlations generated by the dynamics of a system are completely governed by 
cumulants fl20|) of the groups of operators ([7]) of the von Neumann equations. 

The links of constructed solution of the von Neumann hierarchy both with the solution of the 
BBGKY hierarchy and with the nonlinear BBGKY hierarchy for marginal correlation operators 
were discussed. Solution fl^Sjl of the BBGKY hierarchy f lFTI) for Fermi and Bose particles with 
a n-body interaction potential was constructed. This solution is represented in the form of 
expansion over the clusters of particles, which evolution is governed by the corresponding order 
cumulants (13T|) of the groups of operators of the von Neumann equations ([7]). The cumulants 
for the Fermi and Bose systems of particles have the same structure, as in case of the Maxwell- 
Boltzmann statistics [2^] . 

The cumulant structure of solution (122]) of the von Neumann hierarchy (IT^ induces the 
cumulant structure of solution expansion fl62|) both of the initial-value problem of the non- 
linear BBGKY hierarchy for marginal correlation operators and of initial-value problem of 
the BBGKY hierarchy for marginal density operators. Thus, the dynamics of infinite-particle 
systems is generated by the dynamics of correlations. 

As we mention above, for initial data from spaces £^(7"^) the average number of particles is 
finite. In order to describe the evolution of infinitely many particles [21], in particular, to derive 
the quantum kinetic equations from constructed dynamics we have to construct solutions for 
initial marginal density operators or marginal correlation operators that belong to more general 
Banach spaces than Ili(J'^). 
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